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ABSTRACT 

We consider quantum mechanical corrections to a homogeneous, isotropic 
and spatially flat geometry whose scale factor expands classically as a general 
power of the co-moving time. The effects of both gravitons and the scalar 
inflaton are computed at one loop using the manifestly causal formalism of 
Schwinger [|lj with the Feynman rules recently developed by Iliopoulos et al. 
[01. We find no significant effect, in marked contrast with the result obtained 
by Mukhanov et al. 0, f|] for chaotic inflation based on a quadratic potential. 
By applying the canonical technique of Mukhanov et al. to the exponential 
potentials of power law inflation, we show that the two methods produce the 
same results, within the approximations employed, for these backgrounds. 
We therefore conclude that the shape of the inflaton potential can have an 
enormous impact on the one loop back-reaction. 
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1 Introduction 



In recent papers, Mukhanov, Abramo and Brandenberger [Q, |J have inves- 
tigated the problem of the back reaction of quantum fluctuations in chaotic 
inflation. Using canonical quantization, they calculated the one loop effec- 
tive energy-momentum tensor induced by quantum fluctuations of the mat- 
ter and metric fields. They found that the expansion rate of the background 
Friedmann-Robertson- Walker space-time slows down as a result of quantum 
deformations to the equations of motions. 

Is this result generic for all inflationary models? In other words, does 
back reaction slow down the expansion rate H in scenarios other than chaotic 
inflation? In this paper we extend the work of Mukhanov et al. to a very 
different class of inflationary models, the so-called power-law inflation || in 
which the expansion of the universe is given by a scale factor a(t) oc t s with 
s > 1. For that purpose we employ the standard formalism of covariant 
quantization, using the Feynman rules for these backgrounds that have been 
worked out by Iliopoulos, Tomaras, Tsamis and Woodard f2|. 

We find that the one loop back reaction in power-law models is negligible 
and stays negligible through the inflation of the universe. By comparison, 
in the chaotic models the one loop infrared corrections can become non- 
perturbatively large at late times. 

The effect discussed by Mukhanov et al. is closely related to the one 
found in earlier work by Tsamis and Woodard in the context of pure 
gravity with a cosmo logical constant. Both mechanisms have a simple phys- 
ical origin: the enhancement of zero-point energy of the quantum fields by 
the expansion of the universe, also known in the literature as superadiabatic 
amplification. Virtual quanta of cosmological wavelengths become trapped 
by the expansion of the universe and are unable to recombine. Gravitational 
interactions between these virtual pairs, being always attractive, slow the ex- 
pansion rate. Because gravitational interactions are very weak, back reaction 
takes a long time before it can become a sizeable effect. As an infrared effect 
that derives from quantum fluctuations whose wavelengths are comparable 
to the Hubble radius H -1 , back reaction can be studied in the context of 
general relativity, whatever is the ultimate theory of gravityQ 

1 This is precisely the reason why cosmologists are able to predict the seeds of structure 
formation from quantum fluctuations in inflationary models without regard to the true 
theory of quantum gravity^). 
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A second point of this paper is to prove that results obtained by covariant 
quantization should coincide with those worked out using canonical quanti- 
zation. This is relevant to the validity and interpretation of the canonical 
calculation, which was used by Mukhanov et al. to derive their results and 
has recently been challenged by Unruh In a separate paper || we go fur- 
ther and show that covariant quantization yields the same physical results 
as canonical quantization in the case of the chaotic inflation models as well. 

The present paper is organized as follows. In section 2 we describe the 
perturbative background. We write the Feynman rules for the gravity-scalar 
system (taken from Iliopoulos et al 0) in section 3. In section 4 we show how 
to use Schwinger's formalism to obtain expectation values for the metric and 
matter fields from the amputated 1-point functions of the theory. In section 5 
we derive the infrared parts of coincident propagators. Our results are given 
in section 6. In section 7 we re-derive those results using gauge-fixed canonical 
quantization, and show that the two methods give identical results, namely 
that one loop back reaction is negligible in power-law inflation. Section 8 
discusses the implications of this work. 



2 Classical background 

The system we consider is general relativity with a minimally coupled scalar 
field: 

C = T ^RV^9-\d^d^g^^9 ~ V(<p) V=g~ . (1) 
The background is comprised of a flat, homogeneous and isotropic metric 

dsl = —dt 2 + ao 2 (t)dx ■ dx = ^ 2 (rj) (y—drj 2 + dx ■ dxj (2) 

and scalar field <po(t). The conformal time rj is related to comoving time by 
dt = Q(rj)drj. In this paper we examine those scalar potentials which drive 
power-law expansions, that is, the scale factor grows as a function of time 
according to 

ao(t)= (l + f^f , (3) 
where s and H[ are constants. In terms of conformal time, 
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(4) 



and rji = — ^y-H^ 1 . The advantage of this class of expansion laws is that the 
Feynman rules are already known p|. 

We now write some useful formulas for the background quantities. The 
expansion rate is given by the logarithmic derivative of the scale factor with 
respect to comoving time, 



H = !*° - — 
cio Q 2 



(5) 



where a dot indicates d/dt and a prime denotes d/dr] = Qd/dt. 
The Einstein field equations are 



3H 2 
2H - 3H 2 



-K 
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hp\ + V{ip Q ) 
~ V(<po) 



(6) 
(7) 



and the integrability condition for this system is given by the equation of 
motion for the background scalar field, 



(p + 3Hcp + V v (<po) = . (8) 

In the expressions above, k 2 = 167rG is the loop-counting parameter of per- 
turbative quantum gravity, V((p) is one of the scalar field potentials of power- 
law inflation models and V v = dV/dip. 

It is useful to invert the Einstein equations and write the scalar field 
quantities in terms of the Hubble parameter H and its derivative H: 

$ = ~H, (9) 

K z 

V(<p ) = ^(H + 3H 2 ). (10) 
In terms of the conformal scale factor, the expressions are 
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(ii) 

(12) 



If we now substitute the expansion law (U) into expressions flnD and (HD, 
it is easy to obtain the scalar field potentials which correspond to power-law 
inflation ||: 



2 1 

ip = —=-]n(Q) , 
VM = 6(l-l)f exp 



(13) 
(14) 



The parameter s therefore regulates the steepness of the scalar potential as 
well as the rate of expansion. Notice that the equation of state for the energy 
density and pressure of the scalar field is 



V 

w = - 
P 



3H 2 



2 

-1 + 7T- 

3s 



(15) 

-2H-3H 2 3s 
In the limit s —>■ oo we recover exponential inflation and the equation of state 
for de Sitter space, w = — 1. 



3 Quantum theory: Feynman rules 

Before laying down the quantum theory, we would like to address an objec- 
tion that is too often raised: that general relativity is not a perturbatively 
consistent quantum theory of gravity, and so cannot be used to study quan- 
tum effects. 

That judgment is only partially correct. The back reaction of quantum 
fluctuations described here is an infrared effect, caused by fluctuations with 
cosmological wavelengths, and we know, for example from the Fermi theory 
of neutrinos |n|, that infrared physics can be studied by the low-energy ef- 



fective theory, regardless of the renormalizability of that theory. Conversely, 
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whatever the ultraviolet behavior of the true quantum theory of gravity may 
turn out to be, it will not be affected by these infrared effects. As long as we 
are careful not to introduce a spurious time dependence through the ultra- 
violet regularization, back reaction is given at late times by ultraviolet finite 
terms whose form is entirely controlled by the low-energy theory. 

We also note that quantizing some of the gravitational degrees of free- 
dom is crucial to all derivations of the density perturbations and cosmic mi- 
crowave background anisotropies in inflationary cosmology 0. Both we and 
Mukhanov et al. just carry the calculations out to the next stage, and ask 
what effects quantum fluctuations might have on the background in which 
they propagate. 

The goal of this section is to summarize the Feynman rules for the La- 
grangian ([!]). The fundamental fields are the metric and the scalar field, 



g^u = ^ 2 (^ + K^ttO = ^ 2 g^u , (16) 

if = po + <f). (17) 

Here ip^ v is the pseudo-graviton field, whose indices are raised and lowered 
with the Lorentz metric rj^. 

Both the scalar-scalar and the graviton-scalar interaction parts of the La- 
grangian can be easily expanded in perturbation theory from the fundamental 
Lagrangian (|Tj). The pure gravitational interactions are more complicated, 
however, after some integrations by parts the whole Lagrangian reduces to 
the simple expression, 



■ ^wV^] + m'^r a g^ P a,^ou - n 2 ^^g 
i / — / — 00 



1 d n V(<p ) An 



^Ti nl dl P r 



(18) 



up to total derivative terms. 

Gauge fixing is achieved by adding a gauge fixing term and the corre- 
sponding ghost action to the invariant Lagrangian: 
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BRS 



HVTF^F V 



(19) 



The symbol 8F ^ represents the variation of the gauge fixing functional under 
an infinitesimal diffeomorphism parameterized by the ghost field uA We 
follow Iliopoulos et al. in the choice of the gauge fixing functional: 



fl ,1/ 



in' N 



(20) 



A great advantage of this gauge is that it decouples the tensor structure of 
the propagators from their dependence on spacetime. 

The quadratic part of the gauge-fixed Lagrangian (|19"D can be written, 
up to ghost contributions, in terms of a kinetic operator in "super-matrix" 
representation: 



V 



D pa 

fJ,V 



2 i 2s 2 -3s+2 1 



d 2 + 



(s-1) 2 J7 2 



n 



(21) 



where £ M = r^o- To get back the quadratic terms of the Lagrangian (|T9|), just 
multiply the super- matrix on the left by (t/^ 0), and on the right by its 
transpose. The kinetic operator D J* 7 is given by 



D pa 



+ 



D 



B 



(22) 



where barred tensor symbols denotes that the zero components have been 
projected out: 8" = 8" + t^t u . The kinetic operators above are given by 



DS 



n 
n 



d 2 + 
d 2 + 



2s 2 



s 1 



'1 — s) 2 r\ 2 



(1 — s) 2 rf 



n 



n 



(23) 
(24) 



The off-diagonal term in fl2~i"D , coupling ^oo to <f), can be removed by a 
change of variables Q: 
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Cij 

X 

V 



sin 6* + ipoo cos ® 
(j> cos 9 — -000 sin 6 , 



(25) 
(26) 
(27) 



where tan 2 9 = s and Latin letters denote spatial indices. 

In terms of the new fields (25)- (27) the quadratic, gauge-fixed Lagrangian 

is: 



-Si, r<\ 



+ 2~(oi D B(oj [~Sij] 



+ -vD s c v + ^xD s A X + cfDfa + u°D s B u , 



(28) 



where D S A and D% were defined above in (p3|)- (p4|) , and D s c is given by 

2-s 1 



D s c = n 



d 2 + 



1 — s) 2 rf 



(29) 



The quantization of this system is trivial. There are three modes associ- 
ated with the three kinetic operators: 



D' I *(r ] ,k,I) = 



(30) 

We note for now that by Eq. ( p0|) and the definitions of the kinetic operators 
(|23"D -(p4|) and (|29"D , the mode functions are proportional to Hankel functions 
Hj))' 2 \kr]) where the labels vi(s) depends on the parameter s. In the limit 
where kr\ — ► oo we find that the modes proportional to Hff have negative 
frequency, while those proportional to H$ = Hff* have positive frequency. 

The fundamental fields ip an d <p can be expressed in terms of the mode 
functions above. The scalar field, for example, is decomposed in the following 
manner: 



<f>(r),x) 



i 

d 3 k 



(31) 



(2nf 



e l ^(r], k, I)a(k, I) + e"^*^, k, I)a\k, I) 
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The weights fj can be found by inverting relations (|26|)-(|2~T|), so in the ex- 
ample above their values are Ja = sin9, fc = cos 9 and fs = 0. 

The quanta of negative energy are associated with the negative frequency 
modes \l/ and the annihilation operators a(k,I). Therefore the vacuum of 
the theory |0) is defined as the state that is annihilated by all a(/c,/)'s and 
I = A } B,C. 

The creation and annihilation operators a) and a obey the usual commu- 
tation relations 

'a(k,I),a\k',l')] = 5i r 5 {z) (k-k') , (32) 
and the mode functions are normalized by 

#(77, k, I)^'*(r], k, I) - #*(77, k, I)^'(r], k, I) = zfi~ 2 . (33) 
The scalar field propagator is defined by 

tA^(x; x') = (0|T {<j>(x)<P(a/)} |0) . (34) 

We shall need pure graviton and mixed propagators as well. All these can be 
expressed in terms of propagators of the diagonal fields (0\T{ipi(x)ipi(x')}\fy , 
which can be conveniently written in momentum space as 



e - 6 * e *.(*-a) x [9{rj — r] f )^(rj, k, /)#*(?/, k, I) (35) 

+ 6(rf-ri)**(Ti, k, /)*(»/, k,I)} , 

with I = A,B, C . Explicit formulae for these propagators in the infrared 
limit can be found in section 5. 

Finally, we transform from the diagonal variables Qj , v and \ back to the 
original pseudo-graviton ip^ and scalar field cf> by using relations (p5|)-(|27|). 
The set of propagators below is the main result of this section: 



(0\T {ipijixjiprsix')} |0) = iA s A (x;x')2 5i( r 8 s )j - 8ij8 r 

+i [A A (x;x') 

(0\T{M^or(x')}\0) = -iA s B (x;x% 



+i [A A (x; x) + sA s c (x; x')} ——S i:j S rs , 

J. — r S 
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(0|T{^(x)^oo(a;')}|0> 

(0|T{^oo(x)^(^)}|0) 

(O|T{^oo(a:)0(^)}|O) 
(O|T{^(x)0(x')}|O) 

(O|T{0(x)0(x')}|O) 

(0|T{zn (x)5;o(aO}|0) 
(0|T{57 i (x)57 J -(x , )}|0> 



iA^xjx')— + iA^(a;;x') i: ^ , 
(0|T{^oo(^oo(^)}|0)%, 
-iA^x-x^^ + iA^x')^, 
(O|r{^ooW0(^)}|O)^, 

-iA^(x;a/) , 

iA^(x;x')5y • (36) 



4 Attaching external lines 

Back reaction can change the dynamics of the homogeneous and isotropic 
background in three ways: through a redefinition of the time slicing (correc- 
tions to (^oo))) through a change in the scale factor (corrections to the trace 
{ipii)) and through a shift in the scalar field that drives inflation. Because 
the initial state is homogeneous and isotropic, the expectation values of the 
pseudo-graviton and scalar are functions only of time: 



(0|m/v(t7,£)|0> = Ai^ + C^U , (37) 
(0|^,f)|0) = D(rj). (38) 

The expectation values of the full metric and scalar field on the state |0) are 
then 



(0|ds 2 |0) = n 2 [-(1 - C)dr] 2 + (1 + A)dx ■ dx 

= —dt 2 + a 2 (t)dx ■ dx , 
(0M0) = -^=ln(Q)+D, 

V s 



(39) 



(40) 



where t(i]) = J^j. Q(r}')drf is the background comoving time and the scale 
factor in comoving time is given by 
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a 2 (t) = a 2 (t) { 1 + A[rj(t)} + H(t) J dt'C[r](t')] 



(41) 



We don't measure expectation values directly, but rather combinations 
of them that constitute physical observables. The observables should be 
independent of gauge- fixing [[LjJ . In the case of homogeneous and isotropic 
backgrounds, one such physical observable is the effective expansion rate 
H ef[ defined as the logarithmic derivative of the scale factor with respect to 
comoving time: 



H c s = 



d In a(t) 
dt 



H(t) 



(42) 



-c(t) 



H 



A(t) 



2H 2H Jo 



[ dt'CU') + . 
J o 



Rather than computing the 1-point functions A, C and D directly, we 
will compute the amputated 1-point functions instead, and then attach the 
external lines to find the 1-point functions. The amputated 1-point functions 
are defined as 



D pa 



2V? n 2 ±p±c 

s—l rp 

2 _ L 2s 2 -3g+2 1 



d 2 + 



{s-1) 2 rj 2 

<x(jl)y»v + 7(»7) V 
6(r)) 



n 



X 



<0|«^|0) 



(43) 



where the matrix on the left is the same differential operator defined above 
in (PH). Substituting the expectation values Q3"7D and ( P8| ) gives 



5( V ) 



\d A {A - C) , 

Id a (a-c) + d b c-^%d, 



n 



d 2 2s 2 
+ — -r 



3s + 2 1 



dt] 



7} z 



- 1 7] 2 



2^~s tt 2 
s-1?] 2 



c 



(44) 
(45) 
(46) 
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Since we are calculating expectation values, we must employ Schwinger's 
formalism instead of the usual rules for "in-out" matrix elements. For 
the one-loop 1-point functions, the only difference is that in Schwinger's 
formalism the external lines are retarded propagators. In order to enforce this 
choice of external propagator we fix retarded boundary conditions such that 
the 1-point functions A(rj), C{rf) and D(rj), along with their time derivatives, 
vanish on the initial surface rj — rj^. 

The 1-point functions are obtained by attaching the external legs, i.e., 

) with appropriate 



by inverting the coupled differential equations (J44|)-(| 
boundary conditions. After a change of variables and some algebra, one 
obtains 



A 



C 



D 



1 

d2 



i 

d2 



i 

~D~a 



-4a + 
1 



1 (3- + 7)-^' 



D c 
1 

7+1 
l 



s + l 

s 



s + l 



(3a + 7 ) + ^U 



s + l 



(3a + 7 )-^ 



D c 



S (3a + 7 ) + ^U 
s + l 



s + l 



•(3a + 7) + 



s + l 



D 



c 



s + l 



(3a + 7) + 



-5 



s + l 



(47) 



(48) 



(49) 



where Da and Dc have been defined in ( |23| ) and fl29|). Notice that since 
a, 7 and 5 are functions only of time, the non-local operators 1/Dj in the 
expressions above denote two time integrations. The explicit representations 
for the inverted propagators are 



'-J A J Vi J Vi 



(50) 



D A ' 

-±-f(rj) = -n- 2+ -s(r)) fdrf^{rf) f V^ 2+ %W0 , (51) 

L>c Jf]% Jr)i 
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where the lower limits of the integrals make clear that 



A( Vi ) = C( Vi ) = D( Vi ) = , 
A!{jh) = Ci-m) = D'fa) = . 

It is useful to derive explicit formulas for the integrals ([50|)-(|51|) when the 
integrand has the same functional form as the amputated 1-point functions. 
The fastest-growing contributions to the amputated 1-point functions have 
a time dependence oc f2 4_2//s+e (see next section), and for such terms the 
integrals (|5T3l)-(|5T|) reduce to 



lU^) = H^s - 1 + se) (52) 



, r> 3s - 1 + se se _ 3 3 -i 

x -Cl e H Q i 

3s - 1 3s - 1 

„2 



D C K 1 H?(s + 1 + se)(2s - 2 + se 

( ~ f 2s-2 + se„ s+i s + l + se__ 2 ^i 
x -Cl e H Cl 1 



5 Infrared parts of coincident propagators 

In the next section we will sum the diagrams that contribute to the am- 
putated 1-point functions a, 7 and 5. At one loop, the amputated 1-point 
functions are found by acting with the 3-point vertex operators on the propa- 
gators iAj(x; x') and then taking the coincidence limit x — > x'. In this section 
we study the coincidence limits of the propagators and point out those terms 
which can contribute to back reaction. 

The three diagonal propagators can be expressed in terms of the mode 
functions *&(t], k, I) which obey equations (|3(]). After a change of variables 

*(T,,k,l) = {kr,)°h(kr,,I) , w = *L^L , (54) 

we can solve Eqs. fl3"0|) in terms of Hankel functions H Vi '^{krf) of the first 
and second kind, where the indices vi{s) are given in terms of the parameter 
s as 
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VA = \ + JTi ' < 55 > 

, B = i + J i T , (56) 

■*-5-^T- < 87) 

The normalization of the mode functions is given by ([J^). The Hankel 
functions, on their own account, obey the identities 



H^{y) = HW(y), (58) 

HP{y)±H£\v)-H<?\v)±H?{v) = --, (59) 
ay ay 7r y 

which imply that the two linearly independent, normalized mode functions 
are given by 



*(r],k,I) = h\- l ^H^{kr,) , (60) 

<F(T,,k,i) = ±n- 1 ^H!»(h,). (6i) 

The propagators are obtained by substituting the mode functions above 
into the definition (|35f) . In the coincidence limit \\x — x'\\ — > 0, r\ — rf — > 
we have, after performing a trivial angular integration, 

«Aj(x) = dkk*e-* k HM(k V )Hl»(kri) . (62) 

Since H^ ,2 \y) oc y -1 / 2 when y — > oo, the coincident propagators diverge 
as /c 2 in the ultraviolet. After these divergences have been regularized at 
the initial value surface rj = rji, the counterterms should not affect the time 
evolution. Notice that the details of the regularization procedure and the 
ultraviolet behavior of the true, renormalizable quantum theory of gravity 
are issues immaterial to the long-range phenomena described by the effective 
theory, general relativity. 
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We are mainly interested in the infrared behavior of the effective theory, 
since the physical mechanism behind back reaction is superadiabatic ampli- 
fication of quantum fluctuations with physical wavelengths of the order of 
or bigger than the Hubble radius if -1 . We define the infrared as the scales 
larger than the Hubble scale, 

K=^> H(v) = H-.Q 1 ^ . (63) 
In the far infrared the Hankel functions can be approximated by 

(£)" + •••' (64) 

so at coincidence the infrared limit of the propagators (|6^) can be written as 
iA? R \x) = i ^ / dk {knt-*«e-* . (65) 

OTT \T}\\V J 

One immediately sees that by Eqs. Q53|)-(p7D, for large s the integral is 
infrared finite for the B and C modes, but it is divergent for the A mode. 

We can cure this infrared divergence by working on a compact spatial 
manifold. The momenta are then bounded from below, and the integrals 
above should be mode sums. If we set the size of the compact manifold to 
be the H u the Hubble radius at the initial time rji, the infrared cut-off is 
given by ko = H[. In terms of physical wavelengths, the cut-off is the size 
of the initially inflating patch redshifted by the expansion of the universe, 
A, 1 ;'^ 2rr//,0. 

With the cut-off ko = Hi bounding the integral ( p3] ) from below, the 
propagator A is given in the infrared limit (k > HQ) by 

2 

(l - ^- 2/s ) , (66) 

where the numerical factor between square brackets approaches 1 in the limit 
s — ► oo. 

The propagators for the B(+) and C(— ) modes in the infrared limit are 
best left as integrals, 



.c«^_ H l( s -r 



8tt 2 



21 l + ^T fs-l\* 
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4vr 2 



Hi 



(67) 

where it is easy to see that they are dominated by the ultraviolet. In addition, 
they have an overall time factor of Qr 2 l s which makes them subdominant 

( I R) 

when compared with the constant part of zA^ . 

The dominant contribution from the infrared limit of the coincident prop- 
agators is therefore 



2r(| + ^r) ,s-i 



71 



(68) 



6 Results 

In this section we obtain and discuss the results for the expectation values of 
the metric and scalar field in the presence of the quantum fluctuations ip^u 
and 0. 

Our strategy consists of summing up all contributions to the amputated 
1-point functions ([13]) coming from cubic interactions, then obtaining the 
expectation values of the metric and scalar field by d47|)-(|49D. 

There are three types of vertices with cubic interactions in (|T^): pure 
graviton vertices ?/> 3 (Table 1), graviton-ghost vertex ipww (Table 2) and 
vertices with one or more 0's (Table 3). 

For simplicity, we have partially symmetrized the ip 3 vertices of Table 1 
so that the first field always corresponds to the external line. For example, 
the following vertex 

-KH-Stf- 1 ' 3 ^^ , (69) 

gives rise to the symmetrized vertices numbers 1, 2 or 3 of Table 1 when the 
first, third or second ^'s respectively are taken as the external leg. 

As an example, consider the contributions to the amputated 1-point func- 
tion a from the symmetrized vertex number 1. Using the Feynman rules of 
section 3 one obtains 
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M- 
W 


Vei lex raClOI 


M. 
TT 


\ / C 11 ' t OV H O f"f" /ll' 

Vei lex rdLlOI 


i 

1 


±K,Hi\r L / s r) aipi r) a2P2 o 2 t P3> 


22 


1 n2 aofrt^i Bt^Bo o( Q l oft) 
±K,\l 7] 2( 3 7] P3>P2 3 <j{ 


2 


1 r_r o3— 1/s aogo a^iflq j./3i ~l 
^Ktl^r i / s 7 ya2P2^a3P3^ £Pi; 


23 


1 n2 Q-sfai /3i ")(3'< o( Q 2 oft ) 

^KSr?7 a3 ^ Q;i 77 pi ' IP3 aj; o 2 


3 


1 un3-l/s at8i cti B^ o( a 2 j./3o1 

^Kttiir L i s r] azP3 ri \ aipi a\ t P2 > 


24 


1 n2 l( a i SiVa^oft) aoBo 

±K\ro 2 rj Pl>( 3 o 3 i] 212 


4 


—KHilr L / s rj ai ( a2 r] P2 ' pi o 2 ^P3; 


25 


1 n2 o( a 2 fl V m oft) fj-j/3--! 
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Table 1: Vertex factors contracted into V'aiftV'c^ft ^asfo with #1 external. 
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# 


Vertex Factor 


# 


Vertex Factor 


1 


-/tft V 2(Q V l)a3 <9 2 • <9 3 


6 




2 




7 




3 


-Kn 2 rj a2 ( ai d$ l] d? 3 


8 


-2/t J ff i fiV 3( ° 1 '93 l) ^ 2 


4 


2KH i tt 3 r ] a2( - ai d% l) t a3 


9 




5 




10 





Table 2: Vertex factors contracted into ipa^u^u, 



# 


Interaction 


# 


Interaction 


1 




6 




2 




7 


^(3-±)^4-§ #2 


3 




8 


-275( 3 -i)«^ 4_l #^ 


4 




9 


-i(3-i)«fl?n*-^ 


5 




10 


^(3-i)^ 4 -f0 3 






11 


-n=KHjVl? s t n (b „uj^uj a 



Table 3: Cubic interactions involving 0. 
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a {1] (rj) = ixiKX ^KHiQ 3 - 1 / 8 [d a (x)(-iP(x)<ip Qa (x')}] x ,^ x 



2 2s 

-iAa(x; x') -iAc(x; x') 



1 v 7 ' s+1 



(70) 



We are looking for those terms that grow the fastest as a function of time, 
such that after integration by ([52]) or (|53|) there can be a sustained effect of 
back reactions on the metric and scalar field A, C and D. Clearly, the vertex 
above fails to meet that condition and contributes at maximum a constant 
to the expectation values. 

The terms which we are interested in arise when a time derivative hits an 
external leg, as happens in vertex number 3, for example. When that term 
is integrated by parts and the conformal time derivative hits the volume 
factor Hiil 3 ' 1 ^, the result is a factor of if 2 fi 4 ~ 2//s . This sort of term, when 
integrated with respect to the inverse of the A- type propagator by (|52|), gives 
a dominant logarithmic contribution. 

On the other hand, when vertex number 3 is integrated with respect to 
the inverse of the C-type propagator by (|53"D, the result is only a constant 
and other subdominant terms. Since integration by the C-type propagator 
gives only subdominant contributions, we ignore them in what follows. 

It can be verified by inspection that no terms with factors that grow 
faster than f2 4_2 / s arise in the Einstein Lagrangian. Those are the types of 
contributions we are looking for, and only the terms with the right structure 
to contribute a factor of at least if 2 f2 4-2 / s will be collected in the following 
expressions. 

Our results for the amputated 1-point functions to leading order are as 
follows: 



a(rj) = -^(s-i). 2 ^ 4 -/^, (71) 



j{r]) = -a(rj), (72) 
2 



Sfr) = ~a(r,), (73) 



(74) 

where iA^ R ' is given in expression (BH). 
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The expectation values can be found with the help of expressions (f47|)- 
((H) and (||): 



A(rj) 



D(V) 



K 2 Hf s 2 -s 
' 8tt 2 (s + l) 2 

•S 



2r (I + ^ 



S — 1\ 8-1 



X 



lnfi 



7T 

S ^ 3s-l 

+ ~ rfi' 



3s - 1 3s - 1 



— A{ri) 
s 



(75) 

(76) 
(77) 



where we keep the subdominant terms after the logarithm only to stress that 
expectation values vanish at the initial value surface Q = 1. The prefactor 
in square brackets of expression (|75|), which we call <j(s), is approximately 1 
for large s. 

The effective Hubble expansion parameter can be evaluated by substitut- 
ing the expressions above for the expectation value of the metric into Eq. 
([42]). It can be easily seen that the logarithmic contributions cancel, and 
all that remains are subdominant terms. The effective expansion rate is un- 
changed up to terms which are either constant and decaying as functions of 
time, 



H eS (t) = H(t) [1 + (subdominant)] . (78) 

We conclude therefore that there is no one loop back reaction of quantum 
fluctuations on power law inflation to leading order. Presumably there is a 
two loop effect, analogous to the one found by Tsamis and Woodard |J, that 
grows and would become dominant, but we did not attempt to calculate 
those diagrams. 

In the next section we will check the results above with a vastly simpler 
canonical version of this calculation, and provide a physical justification for 
the differences among the power-law and chaotic inflation cases. 
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7 The canonical result 



As a check on our previous results, and in its own right, it is interesting 
to study the same problem of back reaction in power-law inflation using 
canonical quantization. This was the method used by Mukhanov et al. to 
derive their results. 

In order to simplify the full Lagrangian (JI|), we perform the following 
approximations: first, ignore the spin-2 (gravity waves) and spin-1 ("vector" 
fluctuations) projections of the graviton degrees of freedom, and concentrate 
on the scalar fluctuations of the metric that couple to the fluctuations of 
the scalar field (p. Second, use the constraints of the Einstein field equations 
explicitly in (|]) to eliminate the fluctuations of the scalar field and reduce 
the number of degrees of freedom to one: the so-called Newtonian potential 
$. Third, the expectation values are given in terms of the "spectrum" of 
the canonically quantized field $ in power-law inflation, which can be read 
right off the standard formulas in the literature on quantum fluctuations in 
inflationary universe models || [12[ . 



The motivation for the first assumption lies in the observation that the 
infrared limits of coincident propagators of the spin-1 and spin-2 degrees of 
freedom are dominated by the ultraviolet and fall off as a function of time. 
Since cross-correlations are irrelevant at one loop, we can ignore those degrees 
of freedom altogether and concentrate on the dominant, scalar degrees of 
freedom of the graviton. 

This truncated version of canonical quantization has obvious shortcom- 
ings, such as its inadequacy to study perturbative corrections above lowest 
order and the exclusion of many of the degrees of freedom from the calcula- 
tion. These difficulties would be unsurmountable if we wanted to to calculate 
quantum corrections beyond one-loop, for example. Nevertheless, within the 
scope of this one loop calculation the truncated canonical method is perfectly 
suited to study the leading contributions to back reaction. 

The back reaction on the homogeneous and isotropic background is de- 
termined by the expectation values of Einstein's field equations 

(0|GJJ|0> = y(0|ir|0> (79) 
and the equation of motion for the scalar field, 

(o\n<f + vjo) = o, (so) 
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where the inflaton potential is 



V(<p) = 6[l 



3s 





K 




T^P 







(81) 



The terms in Einstein's equations which are quadratic in the quantum 
fluctuations can be collected in an effective stress-energy tensor r Mi ,, which 
is a source term in addition to the stress-energy tensor of the background 
matter || |[. The equation of motion for (p is similarly corrected by source 
terms quadratic in the quantum fields. In this section we shall regard back 
reaction as the response of the background to the source terms induced by 
quantum fluctuations, found by solving for the expectation values of Ein- 
stein's equations and the equation of motion for the scalar field. 

We write the metric in longitudinal gauge, 



ds 2 = - [1 + 2$(t, x)] dt 2 + a 2 (t) [1 - 2$(t, x) + 2w{t)\ dx ■ dx , (82) 
and the scalar field is 

<p = <p (t)+<l>(t,x)+v(t) , (83) 

where $ and d> before, quantum fields which ought to be canonically 

quantized, and w(t) and v(t) are respectively corrections to the scale factor 
and to the scalar field from back reaction^. We have taken advantage of 
the freedom of gauge at the second order to fix any correction to the time 
slicing (the N of earlier sections) to zero. Therefore, in this gauge w is the 
correction to the Hubble expansion rate, H e g = H + w that will eventually 
be compared with the effective expansion rate (f78[ ) of the last section. 

The spectrum \6$(k, t)\ is defined by the expectation value of the square 
of the canonically quantized Newtonian potential: 

(0|$(t,f)$(t,f)|0>= J^\6*(k,t)\ 2 , (84) 

The spectrum is proportional to k when k — > oo, and consequently the 
expectation value above diverges as k 2 in the far ultraviolet. Just as was 

2 We have adopted a notation different from earlier sections to avoid confusion, since 
we in the section we work in a different gauge 
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done in section 5, we assume that a proper renormalization of these infinities 
have been performed and that the appropriate counterterms have no bearing 
on the long-range interactions described by the effective theory. 

The expectation value (|84|) is also divergent in the infrared since the 
spectrum is approximately constant for very small momenta. Again the 
solution is to work in a compact spatial manifold, where k is cut off at the 
value H[ corresponding to the radius of the manifold at the initial time t = 0. 

The infrared limit of the spectrum of the Newtonian potential during 



power-law inflation is [[12 



l4 /R) (M)l 2 



where 



k H? a(s 




(85) 



2r 



i 

s-1 



(86) 



is the same factor defined below (|75|) , and approaches 1 as s > 1. 

The physical information contained in this spectrum is that the ampli- 
tudes of long-wavelength fluctuations are asymptotically constant in power- 
law inflation. The amplitudes of infrared fluctuations in chaotic inflation, 
by comparison, grow slowly with time. Where back reaction is concerned, 
quantum corrections to power-law inflation from a fluctuation with a fixed 
comoving wavelength can be at most constant, while corrections to the back- 
ground in chaotic models can grow as a function of time during the inflation 
of the universe. 

The distinguishing facts about power-law inflation are 1) the equation 
of state w = p/p = — 1 + 2/3s is constant throughout inflation, and 2) the 
kinetic and potential energy densities of the background scalar field, as well 
as the Hubble parameter, are at fixed ratios with respect to each other at 
any given time. Consequently, the amplitude of fluctuations on large scales, 
which couple to these ratios, freeze to a constant value. 

We often find claims in the literature to the effect that quantum fluctua- 
tions during inflation freeze after they become larger than the Hubble radius. 
While this is exact for power-law inflation, it is only approximate for most 
models of inflation. 
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With the spectra above, the expectation value of the Newtonian potential 
in the infrared limit is, by fl34|), 



|$(t, £)$(*, f)|0) = <t(s)'- 



1 k 2 H 2 



32tt 2 



const 



(87) 



Notice that this is in accord with expressions fl36|) and (|68|) for large s (re- 
member that ^oo — 2$). The expectation value of the scalar field can be 
deduced from this expression by using the following useful constraint in mo- 
mentum space, valid in longitudinal gauge: 



<f>{t, k) 



-2^{t,k) 



2^-<£>(t,k) . 



(88) 



The Einstein equations with quantum corrections are found by using the 
metric ( |8"2"D and scalar field ( |8"B| ) into Eqs. ([79]). The general result (the spatial 
gradient terms have been ignored since we consider only infrared modes) is: 



3H 2 + 3($ 2 ) + 12H 2 ($ 2 ) + 6Hw 



K 



2 r 



1 



(89) 



1 

2 



! - ((0 2 ) - 40 o (0$) + 40 2 ($ 2 > + 2<p v) + -V^M) + V„v 



3H 2 + 2H + 4 [H 2 + 2H) ($ 2 ) + 8H($$) + ($ 2 ) + 6Hw + 2w (90) 
[-^o + V - \ ((0 2 ) - 40 o (^) + 40 2 ($ 2 > + 20o*) 



K 

~2 



where brackets (• ■ ■) denote vacuum expectation values. 

As we have argued in the discussion of the spectrum, terms containing 
time derivatives of $ and <fi should vanish. Using the identity (^), we can 
finally write the Einstein equations with quantum corrections above in terms 
of only ($ 2 ), w(t) and v(t): 



6Hw 



K 

~2 



KV 



s ll c-\~[^= + 2(<S> 2 } 



(91) 
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6Hw + 2w 



K 

~2 



KV 



(92) 



It is also instructive to write down the equation of motion for the scalar 
field to second order, 



v + 3Hv + 3wip + V w v (93) 

1, 
2 



- 4(60) - 40 o ($$) + 2V ipip ($<f>) + t,V vw (4> 2 ) = 



which, after use of flS8|), the background identities and the Einstein field 
equation ([H]) read 

t 2 x + -^—tx + 3s (x + 2($ 2 )) = , (94) 

where x = kv / \fs. Since this equation is second order in time, we can always 
find solutions such that v (to) — v{to) = 0. Notice that ($ 2 ) is a constant by 
It is straightforward to solve this equation, and the result is 



KV 

x = —= = —2 



3s-lto 1 fto\ 3s ~ r 



3s - 2 t 3s - 2 V t 



($ 2 ) , (95) 



where we have chosen the integration constants such that at the initial value 
surface x(0) = x(0) = 0. 

The dominant contribution to x — kv/^/s is therefore a constant, cor- 
responding to the homogeneous solution to (p4|). However, as happened in 
the covariant calculation, we have neglected subdominant terms which decay 
as a function of time in the present calculation, by discarding several of the 
degrees of freedom of the fundamental Lagrangian. Therefore, the inhomo- 
geneous solutions are unreliable within our approximation scheme and will 
be discarded. 

Substituting the homogeneous part of (p5|) into (p2|) we get the following 
equation for the correction to the scale factor: 

t 2 w + 3stw = . (96) 
The solution is also straightforward, and we get 
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w = (97) 

exactly since the decaying solution w oc t~ 3s has to be zero to satisfy the 
retarded boundary conditions. 

The effective expansion rate is therefore given by, 

H eS = Hit). (98) 

We have thus obtained the back reaction on the metric and scalar field by 
finding the solutions to the expectation values of the Einstein field equations 
( |79|) and the equation of motion for the scalar field (j80|) . Expectation values 
were evaluated using the canonically quantized Newtonian potential $. 

The result of the canonical calculation is that the back reaction of quan- 
tum fluctuations during power-law inflation does not affect the expansion 
rate of the universe, at least to leading order. This is in agreement with the 
results of the covariant calculation of last section. 

As pointed out earlier in this section, infrared fluctuations of a fixed co- 
moving wavelength in power-law inflation have constant amplitudes. Those 
fluctuations exiting the Hubble radius at a later time in the inflation epoch 
have a smaller amplitude than the ones that exited earlier, since the scale 
of inflation is decreasing as oc t~ 2 . The cumulative effect of superimposing 
modes of different comoving momenta is not sufficient to make the expec- 
tation values of the quantum fields grow in time, and the momentum mode 
sum is dominated by modes that exited the horizon early in the inflation 
epoch. 

8 Discussion 

We have calculated the back reaction of quantum fluctuations on the expan- 
sion rate of homogeneous backgrounds of power-law inflation models. Two 
methods were employed: covariant quantization of the full scalar-graviton 
system, and canonical quantization of a reduced system where the spin-1 
and spin-2 degrees of freedom of the graviton were purged. 

The results of the two calculations are identical: to leading order, there 
is no effect of the quantum fluctuations on the effective Hubble parameter in 
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power-law inflation. The physical reason is that long wavelength modes have 
constant amplitudes, and therefore both the amputated 1-point functions 
and the source terms in the expectation values of Einstein's equations are 
constant. 

The main result is that the shape of the inflaton potential can have an 
enormous impact on back reaction. While chaotic inflation can have a sig- 
nificant back reaction, power- law inflation does not. 

We also hope that the results presented here (see also ||) will help settle 
some of the criticisms raised with respect to the canonical method as applied 
to the study of back reaction in quantum general relativity. 
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